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Abstract 

We define catalytic networks as chemical reaction networks with an essentially catalytic 
reaction pathway: one which is "on" in the presence of certain catalysts and "off" in their 
absence. We show that examples of catalytic networks include synthetic DNA molecular circuits 
that have been shown to perform signal amplification and molecular logic. Recall that a critical 
siphon is a subset of the species in a chemical reaction network whose absence is forward invariant 
and stoichiometrically compatible with a positive point. Our main theorem is that all weakly- 
reversible networks with critical siphons are catalytic. Consequently, we obtain new proofs for 
the persistence of atomic event-systems of Adleman et al, and normal networks of Gnacadja. 
We define autocatalytic networks, and conjecture that a weakly-reversible reaction network has 
critical siphons if and only if it is autocatalytic. 

1 Introduction 

Biological systems exhibit exquisite structure and behavior. We wish to view such sophistication 
in biological systems as the result of, as well as intended for the performance of, computation. It is 
hoped that an attitude to seek for algorithms underlying biological computation will give us insight 
into this sophistication. One of the languages in which the theory of algorithms can be made precise 
is the language of circuits. What are the circuits, if any, that make up biological systems? Over 
the past half century, molecular biologists and biochemists have amassed considerable data about 
biochemical reaction networks. We will postulate that these are the circuits we seek. 

Circuits are constituted from the repeated composition of a small number of distinct parts. 
This circumstance proves of considerable aid in their design and analysis. For example, one may 
view circuits as made up solely of switches. Indeed, the circuits in one of the first computers, the 
Z3, were of this type [5] . What are the switches of biochemical circuits, if such exist? 

To answer this question, let us first consider a catalyzed chemical transformation of species A 
to species B, where species C is the catalyst. In principle, if C is left out, A still gets converted to 
B, but at a lower rate. In practise, the lower rate is typically so much lower that the uncatalyzed 
reaction is left out of the model. This is so especially for biochemical reaction networks — for 
example, one does not expect DNA to self-ligate in the absence of a ligase enzyme within the time 
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scales under consideration, even though such a reaction is possible in principle. We submit that 
catalysts are the switches of biochemical circuits: their presence turns a certain reaction "on," and, 
for all practical purposes, their absence turns it "off." Intuitively, a reaction network is catalytic iff 
there exists an essentially catalytic reaction pathway: one which is "on" in the presence of certain 
catalysts and "off" in their absence. One of the contributions of this paper is to give a mathematical 
formulation of this notion in Definition 13.21 

What dynamical model should we assume for biochemical circuits? There are several choices: 
mass action, chemical master equation, reaction-diffusion, empirical dynamics, etc. Since very little 
is known about the true dynamics in a cell, we should only make the most frugal assumptions. 

The results proved in this paper hold at the reaction network level. Fixing the network con- 
strains, but does not completely specify, the dynamical model. Therefore, our results are robust 
across a wide range of dynamical models, including mass action kinetics. This is in the spirit of 
the approach advocated by Feinberg [TT] . 

Another contribution of this paper is to relate catalysis to the global attractor conjecture, and 
the persistence conjecture, both of which are long-standing open problems in the theory of chemical 
reaction networks. The next few paragraphs will explain this connection. 

In 1974, Horn [T5] made the global attractor conjecture for complex balanced mass action 
systems. Complex balance is a condition that restricts specific rate constants to be special enough 
to guarantee the existence of a "free energy function." Intuitively, the conjecture asserts that 
positive steady states act like global attractors. That is, every solution trajectory which originates 
in the positive orthant must asymptotically reach some positive steady state. This conjecture 
remains open. 

One says that a system of differential equations in Euclidean space is persistent if no solution 
trajectory starting in the positive orthant approaches the boundary. Horn |18j observed that for 
complex balanced systems, persistence and the global attractor conjecture are equivalent. Note 
that there exist mass action systems that violate persistence. For example, the reaction x — ¥ y 
leads to exponential decay of the concentration of species x to zero. 

In 1987, Feinberg |12[ Remark 6.1.E] conjectured that persistence must hold if the underlying 
reaction network is "weakly-reversible": i.e., each component is strongly connected. Since all 
complex balanced systems are weakly-reversible, this conjecture generalizes Horn's conjecture. Very 
little is known about Feinberg's conjecture: it remains open even when all reactions are reversible 
and all reaction rate parameters take the value 1. 

In our hands, the notion of "siphon" turns out to be the connecting link between catalysis and 
persistence. A siphon is a subset of the reacting species whose absence is forward invariant. A 
siphon is critical if and only if the absence of the siphon species is stoichiometrically compatible 
with a positive point. These notions were anticipated by Feinberg |X2(, Proposition 5.3.1, Re- 
mark 6.1.E] without giving them a name. He proved that the zero coordinates of steady state 
points are siphons, and that in the absence of critical siphons, solutions from the positive orthant 
can not asymptotically approach a boundary steady state. Note that this result does not rule out 
trajectories approaching the boundary from the positive orthant. For example, a solution trajectory 
originating in the positive orthant may have an omega-limit on the boundary without violating this 
result. 

The term siphon was known in the literature of Petri net theory, and was introduced to reaction 
networks by Angeli, De Leenheer and Sontag. They proved that the absence of critical siphons 
implies persistence [U Theorem 2]. Their result is non-trivial and requires delicate arguments 
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about omega- limits. 

Anderson and Shiu [21 [3] have analyzed some types of critical siphons for which persistence can 
be established. Shiu and Sturmfels [25J have shown how to compute critical siphons using com- 
puter algebra software. Recently, Gnacadja [H] has analyzed a related condition called "vacuous 
persistence." 

This paper reinterprets critical siphons in the context of biological systems performing com- 
putations. We suggest that to implement exponential amplifiers, single-molecule detectors, and 
molecular logic gates, one requires networks with critical siphons. Our main theorem (Theorem 4.1 ) 
establishes that all weakly-reversible networks with critical siphons are catalytic. One may say that 
the obstruction to proving Feinberg's persistence conjecture comes from catalytic species that act 
like switches. 

To prove the theorem, we show that the network structure that encodes catalysis, and the 
network structure that encodes critical siphons, can both be translated into properties of binomial 
ideals. The crucial insight is that in weakly-reversible networks, the notion of critical siphon does 
not depend on whether individual reactions are reversible or irreversible (Lemma |4.4[ ). This allows 
us to choose all reactions to be reversible, and set all specific rates to be 1. Following [I], this 
detailed balanced mass action system can be represented as a set of binomials. The property we 
seek to investigate manifests itself as a geometric property of the roots of this set of binomials, and 
verifying this completes the proof. 

The proof employs some elementary algebraic geometric ideas (primary decomposition, satu- 
ration). This should not be surprising. Indeed, as has been remarked before [HE], the theory of 
reaction networks is intimately connected with the theory of binomials and binomial ideals. Ex- 
plicitly acknowledging this can sometimes lead to insights, and proof directions, as in the present 
case. 

We introduce the notion of autocatalytic networks (Definition 



6.2), and conjecture that critical 



siphons occur precisely in autocatalytic networks (Conjecture [I). We analyze an example of a 
synthetic DNA molecular circuit, the "seesaw gate" [23], and observe the coincidence of exponential 
amplification, autocatalysis, and critical siphons (Example 4.2). 



2 Preliminaries 

The formal mathematical study of reaction networks was pioneered in the 1970s by Horn, Jackson 
and Feinberg EU EH El EEl ES] . More recent overviews include [fl El E3 ES [26] . 

Intuitively, a reaction network is a graph. Each node represents a complex, or set of chemical 
species, with each species accompanied by a non-negative integer encoding its multiplicity (or 
stoichiometry) . It is conventional in chemistry to represent complexes in additive notation: e.g., 
x + 2y + 5z where x,y,z denote chemical species. Following pQ and Section 1], we will depart 
from this convention and use multiplicative notation for complexes: e.g., the complex x + 2y + 5z 
is now written as xy 2 z 5 . Each edge represents a reaction from the source complex to the sink 
complex. 

The multiplicative notation suggests that reversible reactions can be represented by binomi- 



als p]. This representation is key to the notion of "associated event-system" (Definition 2.6). 

Succinctly, a reaction network is a graph where the nodes are labeled by monomials in the 
species. A reaction network is called weakly-reversible if, whenever there is a path (sequence of 
directed edges) from one complex to another, there is a path back. To assign a dynamics to the 
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network, we need to assign a specific rate constant to each reaction: this is the rate at which 
the reaction would take place if all species had unit concentration. With such an assignment, the 
reaction network is called a mass action system. We recall the definitions, following Section 1]. 

Definition 2.1. 

1. A "chemical reaction network" (CRN) consists of the following data: 

(a) Positive integers s (for number of species) and n (for number of complexes), 

(b) A finite directed graph G with vertices V(G) = {1, 2, • • • , n}, and edges E(G) QVxV, 
and 

(c) An s-variable monomial labeling of the vertices 

(s s s 

11 1 -ll- r f •••• 'II 
j=l j=l j=l 
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such that if ifti = ipy then i = i' ' . 

2. A CRN is "weakly-reversible" iff each connected component is strongly-connected. 

3. A "mass action system" (MAS) is a CRN with a weight function k : E{G) — > M>o- 
Notation 2.2. 

1. The matrix denotes the negative of the Laplacian of G. That is, when i ^ j, the entry in 
position (i,j) is k(i,j), and the row sums are zero. 

2. The matrix of non-negative integers (yij) n xs is denoted by Y = Y^. The i th row (yn,yi2, ■ ■ ■ ,yi s 
is denoted by yi. 

3. The tuple (xi, X2, • • • , x s ) is denoted by x. 

4. For all positive integers m, [m] denotes the set {1, 2, • • • , m}. 

5. The set M s = {J\ s i= iX°[ % \ a £ Z> } denotes all monic monomials in s variables. 

6. The graph Uq is the underlying undirected graph of G with nodes V{G) and edges {(i,j) \ 
either (i,j) E E{G) or (j,i) e E(G)}. 

Definition 2.3 (Mass Action Kinetics). The associated dynamical system of an MAS (G,t/j,k) is 
given by: 

d -^ = ^ x {t)).A k .Y 

We will sometimes write only the graph for CRN's and MAS's. A reference to G could mean 
either the graph G, or the CRN {G,ip), or the MAS (G,ip,k). The meaning should be clear from 
context. 

Definition 2.4. The stoichiometric subspace S = Sg of a CRN G is the linear subspace of M s 
spanned by {yi — y^ \ (i,i') £ E(G)}, the row-differences of Y. 
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For all x, ip(x) ■ Aj~ ■ Y lies within S. Therefore, all solutions to the associated dynamical system 
of an MAS lie within affine translates of S. Moreover, if the solution originates at a non- negative 
point (i.e., x(0) G M> ), then it is known that the solution trajectory is confined to the polyhedron 
obtained by intersecting the appropriate affine translate of the stoichiometric subspace with the 
non-negative orthant [TJ Theorem 4.5], [261 Corollary 7.3]. Following [TJ, this polyhedron will be 
called the "invariant polyhedron." 

Definition 2.5. Let G be a CRN. For all x G IR>oi the invariant polyhedron containing x is 
P x = (s + 5)nM| . 

Note that invariant polyhedra do not necessarily intersect the positive orthant M^q. 

Note that 1 is a monic monomial. Given two monic monomials M = Yii=i X T an d N = Yii=i x { 1 
from M s , recall [I] that M precedes N (and we write M -< N) iff M ^ N and for the least i such 
that ej ^ fa, we have ej < fa. 

Definition 2.6 (Associated Event-System). Let G be a weakly-reversible CRN. The "associated 
event-system" Eg of G is the set of binomials {ipi — ipj \ ipi -< ipj and either belongs to E(JJq) 
or (j,i) belongs to E(U,G)}. 

The monomial order -< chosen in [1] was arbitrary; any total order on monomials would do to 
ensure that the map from a CRN to its associated event-system is well-defined. 

Note that Eg is a "natural event-system" in the language of Adleman et al. [U Definition 2.8]. 
It carries the same information as the MAS {Ug, 1), but this information is expressed in terms of 
polynomials. This is because of the one-to-one correspondence between binomials in Eg and edges 
in U G . 

Recall that a pure difference binomial is a binomial of the form M — N where M, N G Ai s are 
monic monomials. An ideal generated by pure difference binomials will be called a pure difference 
binomial ideal. Note that when G is a weakly-reversible CRN, (Eg) is a pure difference binomial 
ideal. 

We will consider the polynomial ring k[z] = k[xi, X2, ■ ■ ■ ,x s ] where k is a field of characteristic 
zero. (Eg) will denote the ideal generated by Eg- Note that the ideal Zg Q Q[x]/(ei%2 • • ~ x s) of 
Shiu and Sturmfels [53] is the quotient of (Eg) when k = Q. 

Example 2.1. Consider the MAS (G,ip,k): 




The number of species is s = 2. The number of complexes is n = 4. Each complex is labeled 
with a monic monomial in x and y. The CRN is weakly-reversible, since the only component is 



5 



strongly-connected. The associated dynamical system is given by: 



(x,y) = (xy 2 ,x A ,y 3 , 1) 
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(3fei,2^y 2 - 4/c2,33; 4 + ^4,1) -2h,2xy 2 + 3k 2 ,3X 4 - k 3 ^y 3 + 2fe 4i i). 



The stoichiometric subspace Sq is M 2 , and the invariant polyhedron is M 2 . The associated event- 
system Eg is {xy 2 — x 4 , y 3 — x 4 , 1 — y 3 , 1 — xy 2 }. 



3 Catalysis and Catalytic Networks 

The remarkable catalytic behavior of enzymes is one of the most striking features of biochemical 
networks, and arguably central to life itself. Definition |3.2| makes precise what it means for a 



network to behave in a "catalytic" manner. The next three examples motivate Definition 3.2 

For a single reaction, one says that a catalyst is a species whose availability changes the rate of 
a reaction but which is left unchanged by the reaction. 

Example 3.1. For the reversible reaction {x + y „ - x + z}, x catalyzes the conversion of y to z. 

The above example suggests the following rule for identifying catalysts in a system of reactions: 
a species X{ is a catalyst for the CRN G if and only if Xi occurs on both sides of some reaction. 

Example 3.2. Consider the reactions 

x + y ^ ± p 
p + q^=±x + z 



There is no species that occurs on both sides of some reaction. Hence, by the proposed rule, one 
would conclude that there are no catalytic species. However, consider the following sequence of 
reactions: x combines with y to form p, p combines with q to form an x and a z. The net result is 
that x, y and q combine to form an x and a z. The availability of x changes the rate of this reaction 
pathway, but x is left unchanged by the pathway. Hence, it appears desirable to call species x a 
catalyst. 

From Example |3.2[ we see that in the presence of multiple reactions, the concept of catalyst 
requires a little more care in definition. We could postulate that for a network of reactions, a 
catalyst is a species whose availability changes the rate of some reaction pathway, but which is left 
unchanged by that reaction pathway, and attempt to make this precise. Instead, let us turn our 
attention to a related question concerning the entire network. 

It is a common expectation in chemistry that a chemical transformation which takes place in 
the presence of a catalyst must also take place in its absence, though perhaps at a much slower rate. 
However, most biochemical networks of interest violate this expectation. For example, consider the 
ligation of a nicked DNA double strand. In the absence of DNA ligase, this expectation leads us to 
believe that the ligation still takes place. However, in practise, the rate is so low that this reaction 
is never explicitly included in the network. 



6 



Networks that have such "essentially catalytic" pathways we will call "catalytic" networks. 



When framing the formal definition, it is our intention that the networks of Examples 3.1 and 3.2 
be catalytic because, in both cases, there exist reaction pathways on which species x acts like a 
switch. 

Consider a biochemical network that models the ligation of nicked DNA double strands. Suppose 
the very slow reactions corresponding to the spontaneous self-ligation of the DNA double strands in 
the absence of ligase are also included in the network. One expects kinetic simulations to be at least 
as accurate as when these very slow reactions are left out. On the other hand, steady-state analysis 
of such an enhanced network may be much less accurate. This may seem surprising at first, for how 
can a more accurate network model lead to a less accurate prediction? This apparent paradox is 
really an issue of time scales. If it takes an astronomical amount of time to approach steady state, 
then the steady state prediction, though perfectly valid at some time scale, is of doubtful utility. 
Therefore, for the purposes of steady-state analysis, it is something of an art to leave out reactions 
that are not expected to play a substantial role within time scales of interest. The networks that 
omit such very slow reactions are likely to be catalytic, especially in the context of biochemistry, 
because the enzymes that speeds up these very slow reactions acts like switches. 

One more nuance needs to factor into our definition of "catalytic," as we now illustrate. 

Example 3.3. Consider the reactions 

x + y + w ^ ± p + q + w 

We informally identify the species icasa catalyst because its availability changes the rate of the 
reaction x + y + w v ± p + q + w, and w is left unchanged by the reaction. There is no reaction 

of the form x + y ^ - p + q. Therefore, one may be tempted to say that this network is catalytic, 

and that the species w acts as a switch on the reaction pathway from x + y to p + q. However, 
for w to be essential to this pathway, we intend that the absence of w isolate the two ends of 
the pathway. In particular, suppose x and y have positive concentrations, and p and q have zero 
concentration. Then the concentrations of p and q should rise off zero if, and only if, w is present. 
For this network, the parallel pathways x ., - p and y ^ " q together provide a "leakage current" 

which make w non-essential. We will define "catalytic" in such a way as to exclude networks like 
this one. 



Note that in Example 3.3 no vertex in the corresponding CRN is labeled with either of the 
monomials xy or pq. Hence, whether a network is catalytic depends on monomials that may not 
appear as labels of vertices in G. We will extend the graph G to a graph G that includes every monic 
monomial as a vertex, and extends the edges of G in the natural manner. We recall the definition of 
an event-graph from [H Definition 2.9]. For our present purposes, the following simplified definition 
suffices. 

Definition 3.1 (Event-graph). Let G be a weakly-reversible CRN. The event-graph G of G is a 
directed graph with vertices V{G) = Ai s , and edges (Nipi,Nipj) where £ E(G) and N G A4 S 
is a monic monomial. 
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Note that if G is weakly-reversible then every connected component of G is strongly-connected. 
Further, G can be viewed as a subgraph of G in the natural manner. 



Example 3.4. Let G be the network in Example 3.2 The event-graph G has vertices all monomials 



in the variables x,y,p,q, and z. The edges are the obvious edges induced by the network. For 
example, from the monomial xyq, by applying the first reaction, there is an edge to the monomial 
pq. From pq, by applying the second reaction, there is an edge to the monomial xz. 

Definition 3.2 (Catalytic). A weakly-reversible CRN G is catalytic iff there exist path-connected 
M,N G V(G) such that M/gcd(M, A) and N/gcd(M,N) are not path-connected in V{G). 



Example 3.5. Let G be the network in Example |3.2[ The monomials xyq and xz are path- 
connected in G by the path xyq —> pq —> xz. Dividing by gcd(xyq, xz) = x yields the monomials 
yq and z which are not path-connected in G. This last assertion is true because z is an isolated 
node of G, since there is no complex in G that divides z. Therefore, G is catalytic. 

Note that in this definition, gcd(M, A) plays the role of a catalyst. On similar lines, it is easy 



to verify that the CRN in Example |3.1| is catalytic, and that in Example 3.3 is not, as intended. 



Lemma 3.2 will provide an equivalent algebraic characterization of "catalytic" in terms of sat- 
urated ideals. Motivated by this connection, a previous draft of this paper referred to weakly- 
reversible CRN's that are not catalytic as "saturated." 

The content of the next lemma is the link between the ideal {Eg) generated by the associated 
event-system in k[x], and the event-graph G. 

Lemma 3.1. Let G be a weakly-reversible CRN. Let M and iV be distinct monic monomials. Then 
M — N G {Eg) iff M and A are path-connected in G. 

Proof. (<=) Consider a path M = Ml, M 2 , • • • , Mi = N in G from M to A. Then for i = 1 to I - 1, 
Mi - M i+ i belongs to {E G ). Adding, we get M — N £ {£ G ). 

(=>) Suppose M — N £ {Eg)- There exist Cjj& G k and monic monomials Mj,-fc and a positive integer 
K such that: 

M — N = fe c vk {M ijk ■ ip, - M ijk ■ Vi)] • (1) 

(i,j)eE(G) \fc=l / 

Suppose for the sake of contradiction that M, are not path-connected in G. Then M and N 
are vertices in distinct components of G, say K\ and AT 2 . Restricting Equation [l] to monomials in 
V{K\), and noting that for all i,j, k, the monomials M^k ■ ipj and M^k ■ ipi are path-connected in 
G, we have: 

M = ^ c ijk {M ijk • i>j ~ M ijk • ipi) . 

Evaluating this expression at the point (1, 1, • • • , 1), we get 1 = 0, a contradiction. □ 

The next lemma translates the notion of catalytic networks into algebraic terms. Recall from j6l 
Exercise 4.4.8] that the saturation of an ideal I C k[xi,X2, • • • ,x s ] with respect to a polynomial / 
is the ideal I : /°° = {g G k[xi, x 2 , • • • , x s ] \ there exists k > with f k g G /}. 

Lemma 3.2. A weakly-reversible CRN G is catalytic iff {Eg) C {£q) : {x\X2 ■ ■ ■ 
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Proof. Note that the saturation (€g) ■ {x\X2 ■ ■ ■ x s )°° of the binomial ideal (£) is itself a bino- 
mial ideal by [U Corollary 1.7. (a)]. Therefore (Sq) ■ {x\X2 ■ ■ ■ x s )°° is generated by the saturates 
i^mW) I ( M > N ) e E (G)} of binomials. 

(<=) If G is not catalytic then by Lemma 3.1 all binomials ^n^n are in (£g)- 

(=>) Suppose G is catalytic. Then there exist distinct, path-connected Mi,M 2 G V(G) such 

that M 2 /gcd(Mi,M 2 ) to M x j gcd(Mi, M 2 ) are not path-connected in V(G. By Lemma |3~X 



M 2 /gcd(M 1 ,M 2 )-M 1 /gcd(M 1 ,M 2 ) E (£ G ) : ( Xl x 2 - ■ ■ x s )°° \(S G ). □ 



4 Networks with Critical Siphons are Catalytic 

Recall from [1] the concept of "persistence" of positive dynamical systems, which models the idea of 
"species non-extinction." This notion can be extended to chemical reaction networks in a natural 
manner: a CRN is persistent if, for every choice of specific rate constants, the corresponding MAS 
is persistent. Recall that for x : M>o —> M s , the omega-limit set uj(x) of x is the set {y G M s | there 
exists an increasing sequence ti, t 2 , . . . tending to infinity such that the sequence x(t\), x(f 2 ), ... is 
Cauchy with limit y}. In words, the omega- limit set of a trajectory consists of those points that are 
approached arbitrarily close, infinitely often, along some unbounded increasing sequence of times. 

Definition 4.1 (Persistence). A CRN (G,if>) is persistent iff for every MAS (G,ip,k), for every 
solution x : M>o — > K> *° * ne associated dynamical system, the omega-limit set of x does not meet 
the boundary of the positive orthant. That is, u{x) n <9IR>q = 0. 

Note that by our definition a persistent MAS is allowed solutions that escape to infinity, as well 
as "omega-limit points at infinity" where some coordinates are infinite and others may be zero, 
since such points at infinity do not belong to M> , but to a compactification. 

Open 1 (Feinberg's Persistence Conjecture, 1987). \\.2\ Remark 6.1.E] All weakly-reversible CRNs 
are persistent. 

One of the consequences of Feinberg's conjecture is Horn's global attractor conjecture which 
has remained open since 1974. Recall that an MAS (G,ip, k) is complex balanced iff there exists a 
positive point xq G M> such that tp(xo) ■ Af. = 0. 

Open 2 (Horn's Global Attractor Conjecture, 1974). [18] Let G be a complex balanced MAS. 
Then every invariant polyhedron P contains a positive point x* such that for all solutions x(t) to 
the associated dynamical system with x(0) G P n Ml , the limit lim x(t) exists, and equals x* . 

t— >oo 

Recall that an MAS is detailed balanced iff all reactions are reversible, and the specific rate 
constants are such that the reactions admit a positive point of simultaneous balance of all the 
reactions. Every detailed balanced system is also complex balanced. Even for detailed balanced 
systems with every specific rate constant equal to 1, the persistence and global attractor conjectures 
are open problems. 

Recall the notions of siphon and critical siphon from [4J. The intuitive idea is that a siphon 
is a set of species whose absence is forward invariant. A siphon is critical if the absence of the 
siphon species is stoichiometrically compatible with a positive point. The important point to note 
is that this definition can be made at the level of the network, without making any reference to the 
dynamics. 
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Definition 4.2 (Siphon, critical siphon). Let G be a CRN. A nonempty set Z C [s] is a siphon of 
G iff for all (£, j) £ E(G), if there exists k E Z such that x^ \ ipj then there exists I £ Z such that 
xi\ tpi. A siphon Z is critical iff there exists a point z G M> such that Z = {i \ Zi = 0} and the 
invariant polyhedron containing z intersects M^q. 



Example 4.1. For the network of Example 3.1, every siphon containing (the index of) z must 



also contain either x or y. However, the absence of x is forward-invariant, so the set {x} is itself a 
siphon. 

We will now reinterpret critical siphons in the context of biological systems performing com- 
putations. Many circumstances in biology require biochemical circuits to detect extremely small 
concentrations of some species, or the simultaneous presence of several species. This is a signal 
amplification task, and at the same time a logical AND operation. Let us call an operation that 
performs both these tasks simultaneously an "amplifying AND" operation. For example, expo- 
nential amplifiers, single- molecule detectors which must report on the presence or absence of small 
concentrations of some species, and molecular logic gates, all require "amplifying AND" operations. 

For these circuits, the presence of input should lead to amplification of one or more "reporter 
species," preferably exponentially fast so as to exceed a critical threshold of detection within a rea- 
sonable amount of time. We call this the sensitivity requirement for the amplifying AND operation. 
This sensitivity requirement appears to be important for the error-free operation of large circuits. 

In order to ensure that there are no false positives, under appropriate initial conditions, absence 
of input and reporter species must be forward invariant in time. This is a robustness requirement. 
To ensure robustness of the amplifying AND gate, some superset of the reporter species together 
with the inputs (where the extra species act as control) must correspond to a siphon. 

Exponential amplification requires that small and large concentrations of the reporter species 
be stoichiometrically compatible. We suggest that this must lead to a stoichiometrically compatible 
direction transverse to the siphon, and that critical siphons are related to exponential amplification. 
The example below provides evidence in favor of this suggestion. 

Example 4.2. The "seesaw gate" of Qian and Winfree [23] is a synthetic DNA network primitive 
analogous to a transistor. Like transistors, seesaw gates are intended to be composed to form large 
circuits that can perform a variety of functions like amplification, digital logic, etc. Zhang et al. |27| 
have demonstrated exponential amplification with a related motif. 

Here we consider the CRN of a one input, one output seesaw gate, modeled by the two reversible 
reactions below, involving the five chemical species gate : output, input, gate : input, output, and 
gate : fuel. 

gate : output + input ^ ± gate : input + output 
gate : input + fuel ^ - gate : fuel + input 

It is clear that {input, gate : input} is a siphon because the absence of input and gate : input 
is forward invariant. This confers some robustness to a seesaw gate. One wishes to have even 
more robustness against small amounts of noise, since empirical observation suggests that small 
leakages are inevitable when composing larger circuits. Ideally, the output should be released only 
when the input exceeds some threshold. To accomplish this, a small concentration of an additional 
"threshold" species that absorbs "input" species is introduced. This new reaction 

input + threshold ► waste 
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provides negative feedback to small leakages until the threshold species is exhausted. One may say 
that this negative feedback "thickens the siphon." 

Note that the CRN is catalytic. There are edges in the event-graph from the complex fuel + 
gate : output + input to the complex fuel + gate : input + output, to output + gate : fuel + input, 
but there is no path from fuel + gate : output to output + gate : fuel. Thus, input acts as a 
catalyst. 

The siphon {input, gate : input} is not critical because the sum of concentrations of input and 
gate : input is a dynamical invariant, i.e., its value remains unchanged along solution trajectories. 
Correspondingly, one may verify that the seesaw gate is not an exponential amplifier. 

If we set output and input to be the same species, then the net reaction becomes 

fuel + gate : output + input v " gate : fuel + 2 input. 



This leads to exponential amplification, as verified by simulation by Qian and Winfree |23j . Corre- 
spondingly, {input, gate : input} is now a critical siphon. 

For this last reaction network, note that the catalyst input is able to catalyze its own production. 



We will see later (Definition 6.2) that this is an example of an "autocatalytic network." 

We have suggested that critical siphons may be important for achieving interesting behavior in 
biochemical circuits. Our main theorem says that: 

Theorem 4.1. Weakly-reversible networks with critical siphons are catalytic. 



The proof of Theorem 4.1 appears at the end of this section. There are many examples of 
CRNs that have siphons but are non-catalytic. For example, any non-catalytic network that is 
conservative (conserves mass) has a siphon consisting of all the species. Therefore, the adjective 



"critical" is required in Theorem 4.1 



Theorem 4.1 can be used to prove the persistence of non-catalytic, weakly-reversible networks. 
This follows from a non-trivial result of Angeli, De Leenheer, and Sontag, who show that reaction 
networks without critical siphons are persistent [U Theorem 2]. We sketch their proof idea. They 
prove the contrapositive by considering a solution trajectory x(t) that originates in the positive 
orthant (x(0) G M> ) and whose omega- limit set uj intersects the boundary, and produce from uj a 
critical siphon. They consider the set S which is the intersection of this omega-limit set u with the 
boundary <9M> , and claim it is an invariant set. If not forward invariant, then there is a trajectory 
starting on S that enters the positive orthant. Following this trajectory backwards in time, one 
observes that it leaves the non- negative orthant. This is a contradiction, since omega- limit sets are 
invariant. Similarly, we obtain backward invariance. It remains to observe that from this invariant 
set S one can derive a critical siphon, which essentially follows from definitions. 

Though they state their theorem only for conservative networks (where trajectories are bounded), 
their proof nowhere makes use of this assumption. Their result requires neither weak-reversibility 
nor complex balancing. 



Therefore, as a corollary of Theorem 4.1, we have: 



Corollary 4.2. Every weakly-reversible, non-catalytic CRN is persistent. 



In fact, Corollary 4.2 holds for a larger class of dynamical models than mass action kinetics. 



See [1] for details. By Corollary 4.2, weakly-reversible CRNs whose persistence remains unknown 
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must be catalytic. Therefore, one may say that the obstruction to proving Feinberg's persistence 
conjecture comes from catalysis. 



Complex balanced systems are weakly-reversible. Therefore, by Corollary 4.2 non-catalytic 
complex balanced systems are persistent. Since the global attr actor conjecture is true for persistent, 
complex balanced systems [18] , we have: 



Corollary 4.3. The global attractor conjecture is true for non-catalytic, complex balanced MASs. 



We now establish a number of lemmas that will be required for the proof of Theorem 4.1 The 
first task is to translate the notion of critical siphon to algebra. We first translate from the CRN G 
to the underlying undirected CRN Uq which corresponds to setting every reaction to be reversible, 
and from there to the associated event-system £q, which corresponds to setting all reaction rates 
to have the value 1. 

Lemma 4.4. Let G be a weakly-reversible CRN. The CRNs XJq and G have the same siphons, and 
critical siphons. 

Proof. Let ZC [s]. If Z is a siphon of Uq then it is immediate that Z is a siphon of G. 

Suppose Z is a siphon of G. We claim that Z is a siphon of Uq- Consider (i, j) G E(G). We 
need to prove that if there exists k G Z such that \ tpi then there exists I G Z such that xi \ tpj . 
Suppose Xk | tpi- Since E(G) is weakly-reversible, there is a path (j = io,ii, ■ ■ ■ ,i p = i) from j to 
i in E{G). Since Z is a G-siphon, there exists l v -\ G Z such that xi p l \ ipi p _ 1 . By repeating this 
argument backwards along the path till we reach the node j, the claim follows. 

Since G and Uq have the same invariant polyhedra, it is immediate that a siphon Z is critical 
for G iff it is critical for Uq. D 

Note that, by [251 Theorem 3.1], inclusion- minimal siphons of G could have been defined directly 



in terms of the ideal (£g)- This would give an alternate proof of Lemma 4.4 

The critical siphons of Uq can be studied by studying the zeros of £g, as we now show. First 
we need a lemma from linear algebra. 

Lemma 4.5. Let s be a positive integer, S be a vector subspace of M s , b G M^q be a positive point 
and z G S + b be a non- negative point. Consider a G {0, 1} S with coordinates ati = 1 iff Zj > 0, else 
di = 0. Then the invariant polyhedron containing a intersects M> . 

Proof. Let A G M>o be such that for i = 1, 2, . . . , n, \\(bi — z\)\ < 1. Consider d = a + X(b — z). 
Since \{b — z) G S, we have a £ S + d. It is enough to show that d G M^q. For all i = 1, 2, . . . , n, 
either 

Case 1: Qj = 0, in which case Zi = 0, and so di = \{bi — Z{) = Xbi > 0, or 

Case 2: aj = 1, in which case |A(6j — Z{)\ < 1 implies di > 0. □ 

When k = M or k = C, and £ C k[sc], the variety Vk(£) will denote the zeros of £ over k with 
the analytic topology. So far, we have been talking about siphons as subsets of [s]. We will find it 
useful to talk about points whose zero sets are siphons. We call such points siphon-points. 

Definition 4.3 (Siphon-point, critical siphon-point). Let G be a weakly-reversible CRN, and 
Z C [s] be non-empty. A point a G M> is a Z -siphon-point iff a G Vr(£g) and Z = {i \ ai = 0}. 
The Z-siphon point a is a critical Z -siphon-point iff the invariant polyhedron containing a intersects 

W >{y 
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Lemma 4.6. Let G be a weakly-reversible CRN. Z C [s] is a siphon of G iff there exists a Z- 
siphon-point a G M> . Moreover, Z is a critical siphon iff there exists a critical Z-siphon-point 
a. 

Proof. (=^) Suppose Z is a siphon of G. Consider the point a G {0, 1} S with coordinates «j = iff 



i G Z, else = 1. Since Z is a [/^-siphon by Lemma 4.4, it is easily verified that a G Vk(£g)- 



Now suppose Z is critical. From Definition 4.2, there exists a point z G M> such that Z = {i \ 
Zi = 0} and the invariant polyhedron containing z intersects M> . From Lemma 4.5 with S = Sg, 
the invariant polyhedron containing a intersects M>q. 

(<;=) Suppose a G Vr(£q) H M> is such that Z = {i | ai = 0}. Then it is immediate that Z is 
a L^c-siphon. From Lemma |4.4[ Z is a G-siphon. 

Further, if the invariant polyhedron containing a intersects M^q then it is immediate that Z is 
critical. □ 



Our task for proving Theorem 4.1 is now reduced to showing that when a weakly- reversible 
CRN G is not catalytic, Vr{£g) contains no critical siphon-points. By "strongly dense" we will 
mean dense in the analytic topology. 

Lemma 4.7. Let / C C[xi, X2, ■ ■ ■ , x s ] be an ideal such that 1 = 1: (x±X2 ■ ■ ■ x s )°°. Then Vc(I) \ 
Vc(x\X2 ■ ■ ■ x s ) is strongly dense in Vc(I). 

Proof. From [201 Theorem VIII. 3. 5], there exists k G Z>o such that I has a reduced primary 
decomposition I = I± D I2 H • • • Pi I& which is unique up to reordering. In particular, for all j = 1 

to k, h n • • • n n i J+1 n • • • n i k ij. Also, v c (i) = v c (h) u v c (i 2 ) u • • • u v c {h) = 
v c (Vh) u v c (Vh) u • • • U Vc(Vh). 

It is enough to show that for all j, Vc(Ij) \ Vc(xiX2 • • • x s ) is strongly dense in Vc(Ij). Let 
JG{1,2,--- ,k}. 

Suppose Vc(Ij) \ Vc(xiX2 • • • x s ) is empty. Then Vc(Ij) C Vc(xiX2 • ■ ■ x s ). It follows that there 
exists I G Z>o such that {x\X2 ■ ■ - Xs) 1 G Ij. From the reduced primary decomposition of /, there 
exists / such that / ^ Ij and for all i ^ j, f G I{. Then f I but / • (#1X2 ■•■ x s )' G /, a 
contradiction since 1 = 1: (x±X2 • • • x s )°°. 

Hence Vc(Ij) \ Vc(x\X2 • • • x s ) is non-empty. Since Ij is primary, y/Tj is prime and Vc(Ij) is 
irreducible. Further, Vc(Ij) \ Vc(xiX2 ■ ■ -x s ) is an open subvariety of Vc(Ij). The lemma follows 
from [22, Theorem 1, p. 82]. □ 

We proceed by relating VjR(^G)nM> with V<c{£g) and VR(£G)nM> with Vc(£g)\Vc(xiX2 • • • x s ). 

Lemma 4.8. Let G be a weakly- reversible, non-catalytic CRN. Then Vu(£g) H M^q is strongly 
dense in V R (£ G ) n K| . 

Proof. Let p : Vc(£g) ~~ ^ ^>o ^ e given by {zi,Z2," m , z n ) h-> (|^i|, |^|, • • • , |^n|)> Note that the 
image of /> is contained in Vr(£g) because (£g) is a pure difference binomial ideal (i.e., all reactions 
are reversible, and all specific rates are "+1"). Suppose z G Vr(£g) HlR> . We will find a sequence 



in Vm.(£g) H M>q with limit z. From Lemma 3.2, (£g) = (£g) '■ ( x i x 2 ■ • • ^s)°°- From Lemma 4.7 



applied to (£g)j there exists a sequence pi,P2, • • • in Vq(£g) \ Vc(xi%2 ■ ■ ■ x s ) with limit z. Consider 
the sequence p(pi), p(p2), ■ ■ ■ < This sequence lies in R^q because each j>j lies outside Vc(a?iX2 ' ■ • x s ), 
and its limit is z = p(z). □ 
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Proof of Theorem \4-l\ Let G be a weakly-reversible, non-catalytic CRN, let £ = £q, and let 
z £ 9M>q. We will show that z is not a critical siphon-point. From Lemma 4.6 this is enough. 



For all r > 0, let B z {r) = {x G M>q : \\x — z\\ < r} be the intersection with the non- negative 
orthant of the closed ball centered at z and of radius r. For all x G K>Q) let Px denote the invariant 
polyhedron (of Uq) containing x. If P z does not intersect M> > we are done. Suppose P z intersects 
M^q. From continuity, there exists a sufficiently small r > such that for every x G B z (r), P x 
intersects M>q- Let tq be such a sufficiently small r. 

Note that the MAS (Uq, 1) with the weight function k = 1 is detailed balanced: the point 
of detailed balance corresponds to setting each species to have unit concentration. Therefore, for 
every x G B z {tq), from Birch's Theorem [T9], Lemma 4B], the set P x D Vr(£) D M> of detailed 
balanced points in the relative interior of the invariant polyhedron contains exactly one point. For 
every x G B z (ro), let d(x) be the distance of the unique point in P x n Vr(£) H M^q from z. Then 
the function d : B z (ro) — > M>o which sends x to d{x) is a continuous function on a compact set, 
and hence attains its infimum. Let xq G B z (vq) be a point where this infimum is attained. 

If cZ(xo) = then P Xo is not an invariant polyhedron, a contradiction. Hence, d(xo) > 0. Let 
e = min(ro, d g The set B z (e) n Vr(£) n M^q is empty by construction. From Lemma 



4.8 



this 

implies that z £ Vr{£). Hence, by definition, z is not a critical siphon-point. □ 

5 Atoms, Primes, and Catalysis 

The idea that all chemical species are uniquely divisible into "atoms" — species that are immutable 
and indestructible — dates back at least to the work of John Dalton in the early 19th century, 
and perhaps even further back to philosophical speculations in Indian, Greek, and Islamic civi- 
lizations. Mathematical abstractions of this idea in the setting of chemical reaction networks have 
recently been proposed independently by three groups: "atomic event-systems" [T] by Adleman 
et al; "elemented" and "constructive networks" [23j by Shinar, Alon, and Feinberg; and "normal 
networks" and "complete networks" [15J by Gnacadja. Here we present an algebraic formulation 
of these ideas that generalizes all three proposals. In [1], we have proved global stability results 
for detailed balanced atomic event-systems. Gnacadja has also proved global stability for complete 
networks |15j . In this section, we show that our main theorem that critical siphons occur only in 



catalytic networks (Theorem 4.1) generalizes these previous results. Specifically, we obtain persis- 
tence results for atomic event-systems, and normal networks [15] as easy corollaries to our main 
theorem. Constructive networks, on the other hand, may admit critical siphons, and proving they 
are persistent is not amenable to our techniques. Prime ideals will enter into our story, and one 
reason is this easy and well-known algebraic result. 

Lemma 5.1. If I € k[xi, X2, ■ ■ ■ ,x s ] is a prime ideal then for all / ^ I, the ideals 1 = 1: f°°. 

Proof. It is clear that I C. I : f°°. For the sake of contradiction, suppose g G I : f°° \ I. Then there 
exists k > with f k g 6 /. Since I is prime, and g ^ /, it follows that / G /, a contradiction. □ 

Theorem 5.2. Let G be a weakly-reversible CRN. If the ideal [£o] is prime then G is not catalytic. 

Proof. At the point (1, 1, 1), every binomial in £q evaluates to zero, and hence every polynomial 
in (£g) evaluates to zero. But every monic monomial evaluates to 1. In particular, the monic 



monomial £1X2 ■ • • x s does not belong to (£g)- From Lemma 5.1 (£ ) = {£) : (x\X2 ■ ■ • x s )°°. From 



Lemma 3.2 G is not catalytic. □ 
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If / is a binomial prime ideal, then the quotient ring k[xi, x%, ■ ■ ■ ,x s ]/I is isomorphic (though not 
canonically) to a ring of Laurent monomials. In chemical reaction networks, each such isomorphism 
corresponds to a decomposition map of species into "pseudoatoms." We now make this precise. 

The main component of a pure difference binomial ideal I C C[xi, X2, ■ • ■ ,x s ] is the minimal 
prime P containing / such that Vc(I) Pi M> = Vc(P) H M^q. It is well-known that P is also a pure 
difference binomial ideal: it is obtained as the lattice ideal of the saturation of the integer lattice 
corresponding to I [H]. Therefore, k[xi, X2, ■ ■ ■ , x s ]/P is isomorphic to a ring of Laurent monomials. 

Lemma 5.3. Let m,s G Z >0 . For i = 1 to m, let pi = (pa,Pi2,--- ,Pis) and % = (%i, ft.2, • • • ,1is) 
be vectors in Z> . Let £ = {x Pi — x Qi \ i = 1, 2, • • • , m} C C[x]. Let t> i, V2, • • • , vi G Z s be a basis 
for {gi -pi,g2 -P2,-'- ,Qm -Pm} 1 in C s . Let D v : C[x] -> C[ai, ajf 1 , • • • ^a, -1 ] be the ring 
homomorphism that sends X{ h-> a^a^ 1 ■ ■ ■ a? u . Then ker equals the main component of {£■). 

The proof is well-known in the literature, and details can be found, for example, in JSJ. If £ = £-g 
and the decomposition does not use negative numbers (t>i,t>2, • • • ,f/ G Z> ) then G together with 
the choice of isomorphism corresponds to an elemented network. Further, if every "pseudoatom" 
Oj is itself a species — i.e., for each a^, there exists Xj such that Xj t— > Oj — then G is constructive. 
For elemented and constructive networks, the pseudoatoms are referred to as "elements." 

Finite, atomic event-systems were defined in [1] Definition 2.11]. Here we give a chemical 
reaction network version of the definition for convenience, which also generalizes the notion to 
weakly-reversible networks. 

Definition 5.1. Let G be a weakly-reversible CRN. The atoms of G are Aq = {xi \ x% is an isolated 
node in G}. An atomic monomial is a monic monomial all of whose prime factors (variables that 
appear in the monomial) come from Aq. The CRN G is atomic iff every connected component of 
G contains precisely one atomic monomial. 

It is immediate from definitions that pre- normal networks of Gnacadja [15} Definition 6.1] are 
a special case of atomic networks. The definition of pre-complete networks |15l Definition 7.1] 
intuitively says that all species can be decomposed into "elements" (our atoms) in a unique way 
that is preserved by reactions. Further, pre-complete networks are restricted to "reversible binding 
reactions" that are of the form ^ a^Xj ^ - xj, so that every species that is not an element can be 

decomposed further. This gives, for every non-element, a path in the event-graph to a monomial of 
elements. By uniqueness of the decomposition, it follows that pre-complete networks are atomic. 

Lemma 5.4. Let G be a CRN. If G is atomic then (Sq) is prime. 

Proof. Let £ = £q. Let Aq = {01,02, •" i a i} ^ {xi,X2,--- ,x s } be the atoms. We claim that 
k[x]/(£) = k[oi, 0,2, ■ ■ ■ ,aj\. Since the right-hand side is an integral domain, {£ ) must be prime. 

To prove the isomorphism, we give an explicit "decomposition map" D from k[x] to k[oi, <i2, • • • ,ai]: 
Xi gets sent to the unique atomic monomial from Aq that belongs to the connected component 
of Xi in G. By definition of Ag, the atoms get sent to themselves. The map is defined to be the 
identity on k, and extended to polynomials to make it a ring homomorphism. 

If M — N G £ then D(M — N) = 0. This is because M and iV are in the same connected 
component of G, and hence D{M) = D{N). It follows that {£) C kerD. 

To show kerD C (£), first note that kerZ) is generated by pure difference binomials [21, The- 
orem 7.3]. It is sufficient to show that for all pure differences monomial M — N = \\x c ^ — Wx^ 1 , 
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if D(M - N) = then M - N £ (£). We claim that M and D(M) are in the same connected 
component of G. This is certainly true if M = x%, by the definition of D. It now follows in general, 
by the definition of G, and since D is a ring homomorphism. Since D(M) = D(N), it follows that 
the monomials M and ./V are in the same connected component of G. Therefore, by Lemma 3.1 



M-Ne(£). □ 



From Theorem 5.2 atomic and pre-complete networks are non-catalytic, and hence persistent 
due to Corollary 4.2 In this sense, non-catalytic networks generalize the notion of atomicity. 

Recall that every complete network is pre-complete and admits a Lyapunov function. Similarly, 
every finite, natural, atomic event-system is atomic, and admits a Lyapunov function. Therefore, 
we obtain the global attractor conjecture for finite, natural, atomic event-systems [TJ Theorem 6.1], 
and for complete networks \15\ Theorem 8.3], as a corollary to the persistence of non-catalytic 
networks (Corollary |4.2[ ). 

However, our proofs do not extend to constructive networks. Consider the following example. 

Example 5.1. Consider the reaction x ^ - 2x. By choosing the set of pseudoatoms (elements) to 

be the empty set, one gets a map D : C[x] —> C. The map D sends x to 1, and the ring C[x] to 
the ring C. The kernel is (x — 1), which is the main component of the ideal (x — x 2 ). From the 



remark after Lemma 5.3, (G, D) is elemented, with no elements. It trivially becomes a constructive 



network, since there are no elements (pseudoatoms). Correspondingly, {x} is a critical siphon. 



The proof of Lemma 5.4 does not hold for constructive networks because M and D(M) may 
belong to different connected components. In the example above, D(x) is 1, but 1 and x are in 
different connected components of the event-graph. To insist that M and D{M) be in the same 
connected component is to insist that the CRN provide a "decomposition path" from compounds 
to atoms. Constructive networks admit critical siphons because their definition does not insist on 
the existence of such decomposition paths. 

6 Autocatalysis 

The converse of Theorem |4.1| is false. There exist examples of catalytic weakly-reversible CRNs 
with no critical siphons, as we now show. 

Example 6.1. Consider the reversible reactions 

2x^=^0 
x + y^=±y 

There is an edge from the monomial xy to the monomial y. However, there is no path in the 
event-graph from x to 1. Hence, the network is catalytic. However, there are no critical siphons. 
The only siphon is {y} and this is not critical because y is dynamically invariant, i.e., along solution 
trajectories, its value remains unchanged. 

It is possible to make the class of catalytic networks smaller to exclude networks like the one 



in Example 6.1 



Definition 6.1. A weakly-reversible CRN G is strictly catalytic iff there exist path-connected 
M, N £ V(G) such that for all integers k > 0, (M/ gcd(M, N)) k and (N/ gcd(M, N)) k are not 
path-connected in V{G). 
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Note that every catalytic network is strictly catalytic, with k = 1. Weakly-reversible networks 



with critical siphons must be strictly catalytic: Theorem 4.1 Corollary |4.2[ and Corollary 4.3 hold 
with "catalytic" replaced by "strictly catalytic." This is because binomials of the form M k — N k 
can be factorized into the binomial M — N and another factor which is non- vanishing on the positive 
orthant. The non-negative variety Vr(£g) nM> depends only on the binomial M — N. Hence, the 



proofs of Lemma 4.8 and Theorem |4.1| go through even when "catalytic" is replaced by "strictly 
catalytic." 

The next example shows a weakly-reversible CRN that is strictly catalytic but has no critical 
siphons. 

Example 6.2. Consider the reaction x + a ^ - x + b. This reaction is strictly catalytic since 

for all k > 0, there is no path from a k to b k . However, there are no critical siphons. This is 
because x belongs to every siphon, and no siphon-point with value for x is critical, because x 
is dynamically invariant. Adding the reaction x „ - 2x introduces a stoichiometrically compatible 

direction transverse to the siphon, so that x is no longer dynamically invariant, and the set {x} 
becomes a critical siphon. This new network is also strictly catalytic. 

The modification required of the network in Example |6 . 2 1 was of a very special kind. It involved 
the introduction of a reaction where one molecule of x can produce two molecules of x. In other 
words, the species x was acting as an autocatalyst. We now make this precise. 

Definition 6.2. A weakly-reversible CRN G is autocatalytic iff there exist i G [s], a G Z>o, 
M,N G M s such that 

1. xf(M - Xi N) G (Sg) and a* \ M, 

2. For all k G Z >0 , M k - { Xl N) k £ (S G ). 

Every autocatalytic network is catalytic, and in fact strictly catalytic. Note that the last CRN 



of Example 4.2, the CRN of Example 5.1 and the second CRN of Example |6.2| are autocatalytic. 
Here is one more example of an autocatalytic network, where the autocatalysis involves two species. 

Example 6.3. Consider the reactions 

y 2x 
2y^=^x 



The point (0, 0) is a critical siphon-point. The associated event-system is £ = {y — x 2 , y 2 — x}. The 
ideal {£) contains the binomials y — y 4 as well as x — x 4 , however for all k G Z>o, the binomials 
1 — x 3fc and 1 — y 3k do not belong to {£). Hence, the CRN is autocatalytic. 

These examples prompt the following conjecture. 

Conjecture 1. Let G be a weakly-reversible CRN. Then G is autocatalytic iff G admits a critical 
siphon. 

If the conjecture is true, one may hope to attack Feinberg's persistence conjecture (Open[l]) by 
exploiting the special structure of autocatalytic networks. 
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